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Abstract. Astrophysical tests of the stability of fundamental couplings, such as the fine- 
structure constant a, are becoming an increasingly powerful probe of new physics. Here 
we discuss how these measurements, combined with local atomic clock tests and Type la 
supernova and Hubble parameter data, constrain the simplest class of dynamical dark en¬ 
ergy models where the same degree of freedom is assumed to provide both the dark energy 
and (through a dimensionless coupling, to the electromagnetic sector) the a variation. 
Specifically, current data tightly constrains a combination of C, and the present dark energy 
equation of state wq. Moreover, in these models the new degree of freedom inevitably cou¬ 
ples to nucleons (through the a dependence of their masses) and leads to violations of the 
Weak Equivalence Principle. We obtain indirect bounds on the Eotvos parameter r] that are 
typically stronger than the current direct ones. We discuss the model-dependence of our re¬ 
sults and briefly comment on how the forthcoming generation of high-resolution ultra-stable 
spectrographs will enable significantly tighter constraints. 
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1 Introduction 

The discovery of cosmic acceleration at the end of the last century, from luminosity distance 
measurements of Type la supernovas [1, 2], has led to wide-ranging efforts trying to charac¬ 
terize its observational properties and understand its theoretical origin. While a cosmological 
constant A remains the simplest explanation consistent with current observational data, the 
well-known fine-tuning problems associated with this solution imply that alternative scenar¬ 
ios should be sought and actively tested. The most natural alternative explanation would 
involve scalar helds, an example of which is the recently discovered Higgs held [3, 4]. Such 
cosmological scalar helds would lead to dynamical dark energy scenarios. 

If dynamical scalar helds are indeed present, one naturally expects them to couple 
to the rest of the model, unless a yet-unknown symmetry is postulated to suppress these 
couplings [5]. In particular, a coupling of the held to the electromagnetic sector will lead 
to spacetime variations of the hne-structure constant a —see [6, 7] for recent reviews on this 
topic. Indeed there have been some recent indications of such a variation [8], at the relative 
level of variation of a few parts per million, and an ongoing dedicated Large Program at 
ESO’s Very Large Telescope (VLT) is aiming to test them [9, 10]. Regardless of the outcome 
of these studies (ie, whether they provide detections of variations or just null results) the 
measurements have cosmological implications that go beyond the mere fundamental nature of 
the tests themselves. The goal of this work is precisely to address some of these cosmological 
implications. 

Specihcally, in the same spirit of [11, 12], we discuss how astrophysical and local tests 
of the stability of a can be used as additional tests of the underlying dynamical dark energy 
scenarios. We will thus be making the minimal assumption that the same dynamical degree 
of freedom is responsible for the dark energy and the a variations—there are known as 
Class I models in the classihcation of [7]. In this case any observational (astrophysical) or 
experimental (local) tests of the stability of a will directly constrain dark energy. The future 
impact of these methods as a dark energy probe has recently been assessed in some detail 
[13-15], in preparation for forthcoming facilities which include these measurements as a key 
science driver (we will return to these at the end of the paper). Our goal here is to describe 
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how current data already yields useful constraints on Class I models, and thus provide a 
proof of concept for these future facilities. (Any model not in Class I is called a Class II 
model, and we will also briefly comment on them later.) 

In a previous work [16] we obtained a non-trivial constraint on the dimensionless cou¬ 
pling of the scalar field to the electromagnetic sector, at the two-sigma (95.4%) confidence 
level (and marginalizing over the dark energy equation of state wq) 

Id <5x10-6, (1.1) 

significantly improving upon previous constraints, while the ID likelihood for icq (marginal¬ 
izing over () was found to be, at the three-sigma (99.7%) confidence level 

- 1.05 < tco <-0.94. (1.2) 

Both of these were obtained by assuming a fiducial model for the dark energy with a constant 
equation of state, w{z) = wq. Here we will relax this assumption, and study models where 
the dark energy equation of state does vary with redshift. Presumably these will be some¬ 
what more realistic, but we have also aimed to preserve conceptual simplicity by choosing 
parametrizations that do not increase the number of free parameters. Studying and quan¬ 
tifying how the above constraints depend on the choice of models (and priors) is one of the 
goals of the present work. 

Our second goal is to show that, since in these models the new degree of freedom in¬ 
evitably couples to nucleons (through the a dependence of their masses) and thereby leads 
to violations of the Weak Equivalence Principle [17, 18], the new constraints on also corre¬ 
spond to very tight constraints on the Eotvos parameter rj quantifying violations of the Weak 
Equivalence Principle (WEP). Although these constraints are indirect and model-dependent 
(since Class I models are assumed) they are comparable to or even tighter than current direct 
constraints. Thus forthcoming tests of the WEP will also provide important consistency tests 
of these models. 

In what follows we start by briefly reviewing the relation between a varying a and 
dynamical energy in the case of Class I models, as well as the relation to the WEP violations. 
We then list the local, astrophysical and cosmological datasets we will be considering in our 
analysis, and proceed to study two representative classes of models, constraining them with 
the aforementioned data. In particular we discuss how the obtained constraints depend on 
the model and the choice of priors. Einally we present some conclusions, as well as a brief 
discussion of the improvements on the sensitivity of these tests expected from future facilities. 
We also provide a brief discussion, in an appendix, of the usage of the Oklo bound to further 
constrain these models. 

2 Varying a, dark energy and the Weak Equivalence Principle 

Dynamical scalar fields in an effective 4D field theory are naturally expected to couple to the 
rest of the theory, unless a (still unknown) symmetry is postulated to suppress this coupling 
[5, 17, 18]. In what follows we will assume this to be the case for the dynamical degree of 
freedom responsible for the dark energy. Specifically we will assume a coupling between the 
scalar field, denoted (j), and the electromagnetic sector, which stems from a gauge kinetic 
function Bf{(P) 

C^F =. ( 2 . 1 ) 
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One can assume this function to be linear, 

Bpiif) = 1 - M , ( 2 . 2 ) 

(with = SvrG) since, as has been pointed out in [17], the absence of such a term would 
require the presence a (j) ^ —(p symmetry, but such a symmetry must be broken throughout 
most of the cosmological evolution. As is physically clear, the relevant parameter in the 
cosmological evolution is the field displacement relative to its present-day value (in particular 
(pQ could be set to zero). In these models the proton and neutron masses are also expected 
to vary, due to the electromagnetic corrections of their masses; while we will not discuss 
this in detail in the present work, one relevant consequence of this fact is that local tests of 
the Equivalence Principle lead to the conservative constraint on the dimensionless coupling 
parameter (see [6] for an overview) 


IClocall < 10 ^ , 


(2.3) 


while in [11] an independent few-percent constraint on this coupling was obtained using CMB 
and large-scale structure data in combination with direct measurements of the expansion rate 
of the universe. 

We note that there is in principle an additional source term driving the evolution of 
the scalar field, due to a F^B'p term. By comparison to the standard (kinetic and potential 
energy) terms, the contribution of this term is expected to be subdominant, both because 
its average is zero for a radiation fluid and becanse the corresponding term for the baryonic 
density is constrained by the same reasons discussed in the previons paragraph. For these 
reasons, in what follows we neglect this term, which would lead to spatial/environmental 
dependencies. We nevertheless note that this term can play a role in scenarios where the 
dominant standard term is suppressed. 

With these assumptions one can explicitly relate the evolution of a to that of dark 
energy, as in [11] whose derivation we snmmarize here. The evolution of a can be written 


Aa 

a 


a — uq 

«o 


- 1 = C«^(</> - 4>o), 


and defining the fraction of the dark energy density 


(2.4) 


Ptotiz) ~ P4>{z) + Pm{z) ' 


(2.5) 


where in the last step we have neglected the contribution from radiation (since we will be 
interested in low redshifts, z < 5, where it is indeed negligible), the evolution of the putative 
scalar field can be expressed in terms of the dark energy properties and as [19] 




Ml! 


( 2 . 6 ) 


with the prime denoting the derivative with respect to the logarithm of the scale factor. We 
finally obtain 


— iz) = C r Jiilt(z)(l+Wt(z))zpF 

a Jo ^ ^ l + z‘ 


(2.7) 
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The above relation assumes a canonical scalar field, but the argument can be repeated for 
phantom fields [20], leading to 

An' I r] 

— (z) = -C + (2.8) 

a Jo + z 

the change of sign stems from the fact that one expects phantom hied to roll up the potential 
rather than down. 

The realization that varying fundamental couplings induce violations of the universality 
of free fall goes back at least to the work of Dicke—we refer the reader to [21] for a recent 
thorough discussion. In our present context, the key point is that a light scalar held such as we 
are considering inevitably couples to nucleons due to the a dependence of their masses, and 
therefore it mediates an isotope-dependent long-range force. This can be simply quantihed 
through the dimensionless Eotvos parameter rj, which describes the level of violation of the 
Weak Equivalence Principle. One can show that for the class of models we are considering 
the Eotvos parameter and the dimensionless coupling are simply related by [6, 17, 18, 21] 

rj « lO-^C^ ; (2.9) 

therefore, the constraint on obtained in [16] and reproduced in Eq. 1.1 leads to the two- 
sigma indirect bound 

r7<2.5xl0“i^ (2.10) 

which is stronger than the current direct bounds that will be discussed in the following 
section. We emphasize that this relation only applies to Class I models, and as we will see 
in later sections the numerical pre-factor in Eq. 2.10 will be slightly different for different 
models within this class. 

Two interesting examples of Class II models for which the relation between the Eotvos 
parameter and the model’s coupling parameter is different have been described and studied 
in the literature. One is the runaway dilaton class of models of Damour et al. [22], for which 
up-to-date constraints were recently obtained in [23] . Another one is provided by Bekenstein- 
type scenarios such as the so-called Bekenstein-Sandvik-Barrow-Magueijo (BSBM) model 
[24]. Constraints on in this model were recently updated by [25] who found, at the two- 
sigma level and from astrophysical measurements of a alone 

ICbssmI < 2.2 X 10-^ (2.11) 

in this case, using the relation between r] and C originally provided by [24], this constraint 
leads to 

< 5 X 10-^^ (2.12) 

we will briefly return to these models in what follows. In general one can say that in Class 
II models the proportionality factor between and the square of the relevant coupling is 
smaller that 10“^, and therefore the couplings in such models are less constrained by WEP 
tests than those of Class I models. 

3 Available datasets 

Our goal is to constrain dynamical dark energy models coupled to the electromagnetic sector, 
by using the following datasets 
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Redshift Redshift 


Figure 1. Currently available fine-structure constant measurements, with the relative values 
plotted as a function of redshift. The data of [8] is shown on the left panel, while the more recent 
data of Table 1 is shown on the right panel. In both cases the error bars include both statistical 
and systematic uncertainties, added in quadrature. Note the difference in the vertical scales of both 
panels. 

• Cosmological data: we will use the Union2.1 dataset of 580 Type la supernovas [26] 
and the compilation of 28 Hubble parameter measurements from Farooq & Ratra [27]. 
These datasets are insensitive to C- (Strictly speaking, a varying a does affect the 
luminosity of Type la supernovas, but as recently shown in [12] for parts-per-million 
level a variations the effect is too small to have an impact on current datasets, and 
we therefore neglect it in the present analysis.) However, this data does constrain the 
dark energy equation of state wq, effectively providing us with a prior on it. 

• Laboratory data: we will use the atomic clock constraint on the current drift of a of 
Rosenband et al. [28], 

- = (-1.6T2.3) X 10"^^yr“^ (3.1) 

a 

which we can also write in a dimensionless form by dividing by the present-day Hubble 
parameter, 

= (-2.2 ± 3.2) X 10“'^. (3.2) 

Hq a 

This is the tightest available laboratory constraint on a only. Other laboratory con¬ 
straints are weaker and also depend on other couplings too; we leave the discussion of 
these for subsequent work, but the interested reader can find overviews of atomic clock 
tests in [29, 30] 

For the models under consideration this translates into 

1 t 

—— = — S 3f2(^o|l + ^o| ) (3-3) 

Hq a V ^ 

where S denotes the sign of (H-rco), so it is -|-1 for canonical fields and —1 for phantom 
fields. 

• Astrophysical data: we will use both the spectroscopic measurements of a of Webb et 
al. [8] (a large dataset of 293 archival data measurements) and the smaller but more 
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Object 

z 

Aaja (ppm) 

Spectrograph 

Ref. 

3 sources 

1.08 

4.3 ±3.4 

HIRES 

[31] 

HS1549+1919 

1.14 

-7.5 ±5.5 

UVES/HIRES/HDS 

[10] 

HE0515-4414 

1.15 

-0.1 ± 1.8 

UVES 

[32] 

HE0515-4414 

1.15 

0.5 ±2.4 

HARPS/UVES 

[33] 

HS1549+1919 

1.34 

-0.7 ±6.6 

UVES/HIRES/HDS 

[10] 

HEOOOl-2340 

1.58 

-1.5 ±2.6 

UVES 

[34] 

HE1104-1805A 

1.66 

-4.7 ±5.3 

HIRES 

[31] 

HE2217-2818 

1.69 

1.3 ±2.6 

UVES 

[9] 

HS1946+7658 

1.74 

-7.9 ±6.2 

HIRES 

[31] 

HS1549+1919 

1.80 

-6.4 ±7.2 

UVES/HIRES/HDS 

[10] 

QllOl-264 

1.84 

5.7 ±2.7 

UVES 

[32] 


Table 1. Recent dedicated measurements of a. Listed are, respectively, the object along each 
line of sight, the redshift of the measurement, the measurement itself (in parts per million), the 
spectrograph, and the original reference. The first measurement is the weighted average from 8 
absorbers in the redshift range 0.73 < z < 1.53 along the lines of sight of HE1104-1805A, HS1700+6416 
and HS1946+7658, reported in [31] without the values for individual systems. The UVES, HARPS, 
HIRES and HDS spectrographs are respectively in the VLT, ESO 3.6m, Keck and Subaru telescopes. 


recent dataset of 11 dedicated measurements listed in Table 1. The latter include the 
early results of the UVES Large Program for Testing Fundamental Physics [9, 10], 
which is expected to be the one with a better control of possible systematics. Figure 1 
depicts these datasets. 

We use these datasets to constrain the dynamical dark energy models described in the 
following sections. The behavior of a will be determined by Eq.(2.7) for canonical equations 
of state {w{z) > —1) and Eq.(2.8) for phantom equations of state {w{z) < —1). While 
in [16] we assumed a model with a constant equation of state w{z) = wq, here we relax 
this assumption and study models where the equation of state is redshift-dependent, but 
nevertheless characterized by a single parameter wq. This is done in the interest of conceptual 
simplicity and also because the currently available data can only weakly constrain models 
with additional free parameters. In any case, this does not prevent us from studying a fairly 
broad class of models. 

Our main interest is in obtaining constraints on the C“'^^o plane, and for this reason we 
will fix the Hubble parameter to be Hq = 70 km/s/Mpc and the matter density to be Umo = 
0.3 (and further assume a flat universe, so = 0.7). This choice of cosmological parameters 
is fully consistent with the supernova and Hubble parameter data we use. Moreover, we 
have explicitly verified that allowing Hq, or the curvature parameter to vary (within 
observationally reasonable ranges) and marginalizing over them does not significantly change 
our results. (This should be clear from the fact that a parts-per-million variation of a cannot 
dramatically affect these cosmological parameters.) It is clear that the critical cosmological 
parameter here is wq itself, as in Class I models it will be correlated with We therefore 
consider 2D grids of C and wq values, and use standard maximum likelihood techniques to 
compare the models and the data. 

We also list here the available direct constraints on the dimensionless Eotvos parameter, 
quantifying violations to the Weak Equivalence Principle. These stem from torsion balance 
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tests, leading to [35] 


(3.4) 


rj = (-0.7 ± 1.3) X 10"^^ , 
while from lunar laser ranging one finds [36] 

7/ = (-0.8± 1.2) X 10"^^ (3.5) 

Both of these are quoted with their one-sigma uncertainties. 


4 The thawing models of Slepian et al. 


We will start by studying a dark energy parametrization which is relatively simple, in the 
sense that it does not involve any more parameters than we already have. The model was 
recently introduced by Slepian et al. [37]. The Friedmann equation has the following form 


HHz) 


fim(i + z)^ + 


jl + z)^ 
f^m(i + z)^ -|- 


1 + U)Q 


(4.1) 


We will be assuming flat universes, so flm + = Ij and the model is therefore characterized 

by three independent parameters: Hq, (which here we keep fixed as previously justified) 
and wq. Note that the physical interpretation of these parameters is exactly the standard 
one, and in particular, wq is still the value of the dark energy equation of state today. 

The dark energy equation of state has the following behavior 


w{z) 


-1 + (1 + Wo) 


H‘^{z) ■ 


(4.2) 


Note that for high redshifts this always approaches -1, and it diverges from this value as the 
universe evolves, reaching wq today. This is therefore a parametrization for thawing models. 
Apart from its simplicity, this choice of parametrization is also motivated by the recent result 
that if physical priors are used, allowed quintessence models are mostly thawing [38]. Figure 

2 illustrates the behavior of the dark energy equation of state and the fine-structure constant 
for relevant parameter choices. 

We can now proceed to comparing this class of models with the available data. Figure 

3 shows the results of this comparison, separately for the Webb et al. data (top left panel), 
and for the Table 1 data (top right panel)—in both cases, the constraints from the astro- 
physical data are shown by the thin red lines. It’s well known that the Webb et al. data 
is not consistent with the null result, and a weighted mean of the dataset yields a negative 
value, corresponding to a slightly smaller value of a in the past [8]. In our analysis we cor¬ 
respondingly find a one sigma preference for a non-zero coupling (with a negative sign for 
a canonical field, and a positive sign for a phantom field). However, the data is compatible 
with the null result at two sigma. 

On the other hand, the Table 1 data is fully compatible with the null result. It’s 
worth noting that in the former case the reduced chi-square of the best-fit model (for the 
a data alone) is Xmin,Webb = while in the latter case it is Xmin,Table = 

be an indication that some of the uncertainties in the Table 1 measurements have been 
underestimated. For comparison we also show in the bottom panel of Fig. 3, in the same 
scale as before (and also in thin red lines), the local atomic clock constraint of Rosenband et 
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Figure 2. Redshift dependence of relevant parameters in the model of Slepian et al. [37]. Top 
panel: the dark energy equation of state, w{z), for various choices of wq. Bottom left panel: 
Relative variation of the fine-structure constant, Ao/a, for various choices of wq, assuming a coupling 
C = 5 X 10“®. Bottom right panel: Relative variation of the fine-structure constant, for various 
choices of (, assuming wq = —0.95. A flat universe with 17^0 = 0.3 has been assumed throughout. 


al. [28] ; it is clear from the plot that this is currently more constraining than the astrophysical 
measurements. 

In Fig. 3 the cosmological data constraints are shown by the blue vertical lines, while the 
combined (cosmological plus astrophysical, or cosmological plus atomic clock) constraints are 
shown by the thick black lines. The role of the cosmological datasets in constraining wq, and 
thus breaking the degeneracy, is manifest. Naturally, we can obtain tighter constraints 
by combining all the datasets; this is straightforward to do since the Webb et al, Table 1 and 
atomic clock measurements of a are all independent. The results of this analysis are shown 
in Fig. 4. We note that the results of this analysis are fairly similar to those of [16], where a 
model with a constant equation of state, w{z) = wq, was assumed. The main reason for this 
is that the atomic clock bound, which currently dominates the constraints (as can be seen 
from Fig. 3), is only sensitive to the present value of the equation of state (ie, wq itself) and 
not to its evolution. 

Finally, we can also obtain the ID constraint on the coupling by marginalizing over 
the dark energy equation of state wq- The results of this analysis are shown in the left panel 
of Fig. 5. Again we confirm that in the case of the Webb et al. dataset there is a one-sigma 
preference for a non-zero coupling, while in the other cases the null result provides the best 
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Figure 3. One, two and three sigma constraints on the ( — wq plane for the model of Slepian et al. 
[37], from Webb et al. data (top left panel). Table 1 data (top right panel) and the atomic clock bound 
(bottom panel). In each panel the thin red lines correspond to the constraints from the astrophysical 
or clock data alone, the blue vertical ones correspond to the cosmological data (which constrain wq 
but are insensitive to C) and the black thick lines correspond to the combined datasets. 


fit. The full dataset allows us to obtain a non-trivial constraint on At the two-sigma 
(95.4%) confidence level we find 


iCsGzl <5.6 X 10-6, (4.3) 

which leads to a constraint on WEP violations 

ilSGZ < 3.1 X 10“^'^. (4.4) 

These constraints are very slightly weaker than those obtained in [16] for the constant equa¬ 
tion of state model (cf. Eqs. 1.1 and 2.10 respectively). Physically, the reason for this is that 
in a thawing model with a given wq the amount of a variation at a given non-zero redshift 
will be slightly smaller than that in a constant equation of state model with the same wq. In 
any case, our indirect WEP bound is still stronger than the available direct bounds, cf. Eqs. 
3.4-3.5. 
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Figure 4. One, two and three sigma constraints on the C. — wq plane for the model of Slepian et 
al. [37], from the full dataset considered in our analysis: Webb et al. data plus Table 1 data plus 
atomic clock bound plus cosmological (Type la supernova and Hubble parameter) data. The reduced 
chi-square of the best fit is Xmin full — 0-97. 



Figure 5. ID likelihood for ^ marginalizing over wq (left panel) and for wq marginalizing over C (right 
panel), for the model of Slepian et al. [37]. Plotted is the value of for cosmological 

-I- Webb data (blue dashed), cosmological -I- Table 1 data (blue dash-dotted), cosmological -|- atomic 
clock data (red dotted) and the combination of all datasets (black solid). 


We can similarly obtain the ID likelihood for wq by marginalizing over (^; this is shown 
on the right panel of Fig. 5. In this case we hnd at the three-sigma (99.7%) conhdence level 

- 1.05 < tco <-0.92, (4.5) 

which is again slightly weaker than the one for the constant equation of state model (cf. Eq. 
1.2). Note that, as can be seen in the blue contours in Fig. 3, from our cosmology data 
alone we would get at three-sigma —1.11 < wq < —0.87, so the improvement provided by 
the OL data is signihcant. Therefore this is nominally a very strong bound, though we note 
that it should be interpreted cautiously, both due to our assumptions on other cosmological 
parameters and also because the likelihood is clearly not Gaussian near the minimum. This 
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raises the issue of the choice of priors—specifically, of the choice of a flat prior for 1 + rco- 
We will discuss this point in the following section. 


5 A class of freezing models 

In the previous section we provided constraints on thawing models. Here we will consider the 
opposite scenario: that of freezing models where the dark energy equation of state evolves 
towards —1. Our motivation here stems from the fact that in many dilaton-type models the 
scalar field depends logarithmically on the scale factor 

(/>(z) oc log (1 + z). (5.1) 


Therefore, for a linear gauge kinetic function as we are assuming here, it follows that in that 
case ^ oc ln(l + z). We will presently calculate the condition on the dark energy equation 
of state for Class I models to have a such a behavior for a(z), but it is worth emphasizing at 
this point that some Class II models are also known to display such a behavior. Indeed, for 
runaway dilaton scenarios this behavior is even approximately true after the onset of dark 
energy domination (see for example [22, 23]). By contrast, in BSBM models (another Class 
II model) the field departs from this behavior and freezes quite abruptly at this epoch [24]. 

One may therefore ask what kind of dark energy equation of state would lead to this 
behavior. From Eq. 2.7 we infer that the function inside the square root therein must be a 
constant, that is 

Vl(j,{z)[l + w{z)] = const .; (5-2) 

this can be recast into the following equation 


dw 

dz 


-3(1 


1 -k rc 
1 -k tco 



(5.3) 


Note that the initial condition for the first derivative is 


dw 

dz 


z=0 


-3^rnWo{l + Wo), 


(5.4) 


and for the second one we could also write 


= SQrnWoil + u;o)[l -k 3wo + 311^(1 + W^o)], (5.5) 

2=0 

so w' ~ 311^(1 + wq) and w" ~ 614^(1 -k wq) near the ACDM limit. 

The above equation can be easily integrated, leading to the solution 



w{z) = 


[1 - 0^(1 -k Wo)]wo 


12^(1 -k rCo)(l + 2)^[^ 


(5.6) 


where as usual we are assuming that 12^ + = 1. An analogous solution was obtained, in 

a different context, in [19]. 

The Friedmann equation in this case has the explicit form 


hHz) 

Hi 


— f2m(l-k2)^-k 




0 


Hm(l -k Wo) - Wq 


11^(1 + u;o)(l + ^)^+ , (5.7) 
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Figure 6. Redshift dependence of relevant parameters in the dilaton-like class of models. Top 
panel: the dark energy equation of state, w{z), for various choices of wq. Bottom left panel: 
Relative variation of the fine-structure constant, Ao/a, for various choices of wq, assuming a coupling 
C = 5 X 10“®. Bottom right panel: Relative variation of the fine-structure constant, for various 
choices of assuming wq = —0.95. A flat universe with Umo = 0.3 has been assumed throughout. 


and naturally the evolution of a is given by 

Arv / 

— (z) = C y 30<^(1-ktuo) ln(l-h^). (5.8) 

Figure 6 shows the behavior of the dark energy equation of state and the fine-structure 
constant in this class of models for relevant parameter choices. 

We could treat this parametrization phenomenologically, allow for values of tco < — 1 
(with a flat prior on 1 -|- wq), and fit it to our datasets as was done in the previous section. 
This would lead to slightly tighter constraints, for the reason already explained: in a freezing 
model with a given wq the amount of a variation at a given non-zero redshift will be slightly 
larger than that in a constant equation of state model with the same wq- The full dataset 
constraints on the C ~ '>^0 plane are shown in Fig. 7, which should be compared to the 
analogous plot for thawing models, Fig. 4. As for ID marginalized likelihoods, these would 
now become |C| < 4.6 x 10“® at the two-sigma confidence level and —1.04 < wq < —0.96 at 
the three sigma conhdence level. 

However, the above analysis may be too simplistic, since in this case one physically 
expects that wq > —1. We therefore discard the phantom part of this parameter, and use 
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Figure 7. One, two and three sigma constraints on the C~'^o plane for dilaton-like models, from the 
full dataset considered in our analysis: Webb et al. data plus Table 1 data plus atomic clock bound 
plus cosmological (Type la supernova and Hubble parameter) data. The reduced chi-square of the 
best fit is Xmin full — 0-97. For this analysis we phenomenologically allowed for values of wq < —1. 
Compare to the analogous plot for thawing models. Fig. 4. 


this model as a testbed for the effects of the choice of priors; instead of the flat prior on 
1 -b tco we have used up to this point, we will now assume a logarithmic one. 

Figure 8 shows the results of this analysis for the individual a datasets, alone and in 
combination of the cosmological data. Again we see that for the Webb et al. data there 
is a one-sigma preference for a non-zero value of the coupling while the data of Table 1 is 
compatible with the null result. The reduced chi-square values for the best-fit model, for the 
a data alone, are now xlrin,Webb = 1-05 and xLn, Table = 1-28 respectively. 

These plots, as well as the corresponding constraints from the combined datasets shown 
in Fig. 9, also make it clear that for a value of wq sufficiently close to rco = — 1 any value of 
the coupling would in principle be allowed—although in practice the local WEP constraints 
should of course be satisfied. In principle, and given the form of Eqs. 2.7-2.8, exactly the 
same would also happen in the orthogonal direction (for a sufficiently small C any wq would 
be allowed), but in practice this is prevented by the strong priors on wq coming from the 
cosmological datasets. 

The ID marginalized likelihoods for and wq are depicted in Eig. 10. The bounds on 
the dark energy equation of state are unaffected, as we again find 



Wq < —0.96 

(5.9) 

at the three sigma confidence 
alone, cf. Eig 8), while those 
confidence level we now have 

level (as compared to wq < —0.92 from the cosmology data 
on the coupling are somewhat weakened. At the one-sigma 


Cdil < 6 X 10 ® , 

(5.10) 

while at the two-sigma level 

ICdilI < 2.5 X 10-^ 

(5.11) 
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Figure 8. One, two and three sigma constraints on the C — wq plane for dilaton-like models with 
a logarithmic prior on the dark energy equation of state, from Webb et al. data (top left panel), 
Table 1 data (top right panel) and the atomic clock bound (bottom panel). In each panel the thin 
red lines correspond to the constraints from the astrophysical or clock data alone, the blue vertical 
ones correspond to the cosmological data (which constrain wq but are insensitive to ^) and the black 
thick lines correspond to the combined datasets. 


translating these into WEP bounds, our one sigma constraint is still stronger than the direct 
bounds, cf. Eqs. 3.4-3.5: we have at the one-sigma conhdence level 

IryoiLl <4x 10-1^ (5.12) 

We thus conclude that although our constraints do exhibit some model dependence (both in 
terms of the class of models being assumed and in terms of the underlying priors), they are 
generically competitive with other existing tests of these models. 

As previously mentioned, in this class of models the logarithmic dependence of the 
scalar field, and therefore that of a, is assumed to persist until the present day. Implicitly 
this assumes that the held is not signihcantly slowed down by the onset of dark energy 
domination. It is instructive to compare our constraints with those obtained in models 
where the held is signihcantly slowed down at this epoch. The BSBM class of models is 
the simplest toy model where this occurs. It can be shown that in this case the drift of the 
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Figure 9. One, two and three sigma constraints on the ( — wq plane for dilaton-like models with a 
logarithmic prior on the dark energy equation of state, from the full dataset considered in our analysis: 
Webb et al. data plus Table 1 data plus atomic clock bound plus cosmological (Type la supernova 
and Hubble parameter) data. The reduced chi-square of the best fit is = 0.97. 



Figure 10. ID likelihood for ^ marginalizing over wq (left panel) and for wq marginalizing over C, 
(right panel), for dilaton-like models with a logarithmic prior on the dark energy equation of state. 
Plotted is the value of ~Xmim cosmological -I- Webb data (blue dashed), cosmological -|- 

Table 1 data (blue dash-dotted), cosmological -I- atomic clock data (red dotted) and the combination 
of all datasets (black solid). 


fine-structure constant at the present day is to a good approximation given by [39] 


1 a 

Hq a 


ZTT 


(5.13) 


Unlike Eq. 3.3 which applies to the models we have been studying, this one does not depend 
on the dark energy equation of state, the reason being that in these models the a variation 
is crafted onto the theory in such a way that it does not significantly affect the background 
dynamics (and therefore a cosmological constant is still assumed to provide the dark energy). 
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In this case the Rosenband et al. bound given by Eq. 3.2 leads to 

Cbsbm = (0.8 ±1.2) X 10-4 (5.14) 

at the one-sigma confidence level. As expected, for the BSBM scenario the fact that the held 
dynamics is strongly damped close to the present time implies that the atomic clocks bound 
is much weaker than the one obtained in [25] from the astrophysical measurements of a, cf. 
Eq. 2.11. 

6 Conclusions and outlook 

In this work we used a combination of astrophysical spectroscopy and local laboratory tests 
of the stability of the hne-structure constant a, complemented by background cosmological 
datesets, to constrain the simplest examples of Class I dynamical dark energy models (in the 
terminology of [7]), where the same degree of freedom is responsible for both the dark energy 
and a variation of a. In these models the behavior of a depends both on a fundamental physics 
parameter (the dimensionless coupling C, of the scalar held to the electromagnetic sector) and 
on the background dark cosmology parameters. For the classes of models we have studied 
these are the dimensionless dark energy density and the dark energy equation of state 
wq, the latter being the crucial one. 

We have obtained new, tighter constraints on the dimensionless coupling and also 
studied their model dependence. Our results indicate that current constraints on Class I 
models will be dominated by the atomic clock bound of [28]. However, this is not the case 
for Class II models—the BSBM model provides an explicit example of this point. The 
constraints are somewhat dependent on the choice of (hat or logarithmic) prior for the dark 
energy equation of state, but regardless of this choice these constraints are non-trivial and 
competitive. 

We note that different currently available astrophysical measurements of a —specihcally 
the archival data of Webb et al. and the dedicated measurements of Table 1— lead to some¬ 
what different constraints, with the former leading to a mild (statistically not signihcant) 
preference for non-zero couplings. For the combined measurements we always hnd results 
consistent with the standard paradigm. In any case these discrepancies highlight the impor¬ 
tance of obtaining improved astrophysical measurements of a (both in terms of statistical 
uncertainty and in terms of control over possible systematics), not only for their own sake but 
also because there can have a strong impact on dark energy studies. The ongoing UVES Large 
Program should further improve the status quo, and the next generation of high-resolution 
ultra-stable spectrographs such as ESPRESSO (due for commissioning in late 2016) and 
ELT-HIRES will be ideal for this task. A broad roadmap for these studies is outlined in 
[7], and some specific forecasts of the future impact of these measurements may be found in 
[14, 15]. 

In the classes of models under consideration the new degree of freedom inevitably couples 
to nucleons (through the a dependence of their masses) and leads to violations of the Weak 
Equivalence Principle. We have therefore used our bounds on C, to derive indirect bounds 
on the Eotvos parameter r]. Despite the aforementioned model dependence, our indirect 
bounds are stronger than the current direct ones, in some cases by as much as one order of 
magnitude: in other words, they are at the rj ~fewx IQ-^^ level. We note that the forthcoming 
MICROSCOPE mission, currently scheduled for launch in April 2016, should reach r] ~ IQ-^^ 
sensitivity [40]. Should this measure a value of rj larger than that in our bounds, this would 
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rule out the Class I models we have studied here (or alternatively would imply that the 
measurements of a on which they rely are incorrect and dominated by systematics). 

Finally, we note that the forthcoming generation of high-resolution ultra-stable spec¬ 
trographs will also provide significantly tighter constraints on jy. Specifically, for Class I 
models and based on the forecasts of [14], we may expect a sensitivity of r\ ~fewxlO“^® for 
ESPRESSO and ry ~ 10“^® for ELT-HIRES. The latter is similar to the expected sensitivity 
of the proposed STEP satellite [41]. Therefore null results from STEP and the E-ELT would 
force any putative coupling of light scalar fields to the standard model to be unnaturally 
small, implying that either WEP violating helds do not exist at all or that these couplings 
are suppressed by some currently unknown symmetry mechanism —whose existence would 
be as exciting and significant as that of a light scalar field itself. 

A Appendix: Including Oklo 

The Oklo natural nuclear reactor also provides a complementary probe of the stability of 
fundamental couplings. In particular, it nominally provides a strong constraint on a, but 
it only does so if one assumes that everything else is not varying. This is very a poor 
assumption, as has been amply documented in the recent literature. We refer the interested 
reader to a recent review on the subject [42] and references therein. It is clear that this is 
not as ’clean’ and reliable a measurement as the atomic clock and QSO measurements, and 
for this reason we have not used it in the main text of this work. 

Nevertheless, one may legitimately ask what impact the a constraint on Oklo, taken at 
face value, would have on our analysis. This constraint, from the analysis of [43], is 

— = (0.5±6.1) X 10-®, (A.l) 

a 

at an effective redshift zokio = 0.14. This nominally strong bound ultimately exploits the 
presence of a 97.3 meV resonance in the neutron capture by the Samarium-149 isotope 
(whereas the typical energy scale of nuclear reactions is of order MeV). We note that even 
stronger bounds have been obtained in [44, 45], but these rely on additional assumptions; in 
what follows we will use the more conservative bound of [43]. 

That said we can revisit the analysis done in Sect. 4 for the model of Slepian et al, 
now also including the Oklo constraint. The thawing class of models (of which this model 
is an example) is the one for which Oklo should have a larger relative effect, since in these 
models the deviations from w = —1 are larger at low redshifts. In freezing models (or even 
the model with a constant equation of state {w{z) = wq) the impact of Oklo will be smaller. 

Figures 11 and 12 depict the Oklo impact: we show the 2D and ID marginalized like¬ 
lihoods for ( and wq, with and without adding Oklo to our datasets. Table 2 compares the 
derived constraints in the two cases. One sees that the effects of the Oklo constraint are 
certainly discernible, but by no means dramatic. The qualitative reason for this is that this 
is a constraint at very low redshift (as compared to the astrophysical measurements), and 
only a factor of 3 stronger than the atomic clock constraint. Thus a reasonable fraction of 
models that £t the latter constraint also fit the former. 
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Figure 11. One, two and three sigma constraints on the ( — wq plane for the model of Slepian et 
al. [37], for the full dataset considered in the main part of our analysis (solid black lines, cf. Fig. 4) 
and for the same dataset augmented by adding the Oklo bound (red dashed lines). In both cases the 
reduced chi-square of the best fit is Xmin full — 0-97. 



Figure 12. ID likelihood for ^ marginalizing over wq (left panel) and for wq marginalizing over C 
(right panel), for the model of Slepian et al. [37]. Plotted is the value of for the 

full dataset considered in the main part of our analysis (solid black lines, cf. Fig. 5) and for the same 
dataset augmented by adding the Oklo bound (red dashed lines). 


Parameter 

Confidence level 

Without Oklo 

With Oklo 

Coupling 

95.4% 

ICsGzl < 5.6 X 10-*^ 

ICsGzl < 4.5 X 10-*^ 

Eotvos 

95.4% 

r/sGZ < 3.1 X 10-14 

VSGZ < 2.0 X 10“^^ 

Eq. of State 

99.7% 

—1.05 < wq < —0.92 

— 1.04 < Wq < —0.93 


Table 2. Constraints on the relevant parameters of the model of Slepian et al. [37]. The middle 
column shows the constraints discussed in Sect. 4, which do not include the Oklo bound. The 
right column shows the same constraints also including the Oklo constraint of [43], interpereted as a 
constraint on a. 


- 18 - 





















PTDC/FIS/111725/2009 (FCT, Portugal). CJM is also supported by an FCT Research 
Professorship, contract reference IF/00064/2012, fnnded by FCT/MCTES (Portugal) and 
POPH/FSE (EC). M.P. and M.v.W. acknowledge financial support from Programa Joves i 
Ciencia, funded by Fundacio Catalnnya-La Pedrera. 

References 

[1] A. G. Riess et al. [Supernova Search Team Collaboration], Observational evidence from 
supernovae for an accelerating universe and a cosmological constant, Astron. J. 116 (1998) 

1009 [astro-ph/9805201]. 

[2] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Measurements of Omega and 
Lambda from 42 high redshift supernovae, Astrophys. J. 517 (1999) 565 [astro-ph/9812133]. 

[3] G. Aad et al. [ATLAS Collaboration], Observation of a new particle in the search for the 
Standard Model Higgs boson with the ATLAS detector at the LHC, Phys. Lett. B 716 (2012) 

1 [arXiv:1207.7214 [hep-ex]]. 

[4] S. Chatrchyan et al. [CMS Collaboration], Observation of a new boson at a mass of 125 GeV 
with the CMS experiment at the LHC, Phys. Lett. B 716 (2012) 30 [arXiv:1207.7235 [hep-ex]]. 

[5] S. M. Carroll, Quintessence and the rest of the world, Phys. Rev. Lett. 81 (1998) 3067 
[astro-ph/9806099]. 

[6] J. P. Uzan, Varying Constants, Gravitation and Cosmology, Living Rev. Rel. 14 (2011) 2 
[arXiv:1009.5514 [astro-ph.CO]]. 

[7] C. J. A. P. Martins, Fundamental cosmology in the E-ELT era: The status and future role of 
tests of fundamental coupling stability, Gen. Rel. Grav. 47 (2015) 1, 1843 [arXiv:1412.0108 
[astro-ph.CO]]. 

[8] J. K. Webb, J. A. King, M. T. Murphy, V. V. Flambaum, R. F. Carswell and 

M. B. Bainbridge, Indications of a spatial variation of the fine structure constant, Phys. Rev. 
Lett. 107 (2011) 191101 [arXiv:1008.3907 [astro-ph.CO]]. 

[9] P. Molaro, M. Centurion, J. B. Whitmore, T. M. Evans, M. T. Murphy, 1. 1. Agafonova, 

P. Bonifacio and S. D’Odorico et al., The LIVES Large Program for Testing Fundamental 
Physics: I Bounds on a change in alpha towards quasar HE 2217-2818, Astron. Astrophys. 555 
(2013) A68 [arXiv:1305.1884 [astro-ph.CO]]. 

[10] T. M. Evans, M. T. Murphy, J. B. Whitmore, T. Misawa, M. Centurion, S. D’Odorico, 

S. Lopez and C. J. A. P. Martins et al., The UVES Large Program for testing fundamental 
physics HI. Constraints on the fine-structure constant from three telescopes, Mon. Not. Roy. 
Astron. Soc. 445 (2014) 1, 128 [arXiv: 1409.1923 [astro-ph.CO]]. 

[11] E. Calabrese, E. Menegoni, C. J. A. P. Martins, A. Melchiorri and G. Rocha, Constraining 
Variations in the Fine Structure Constant in the presence of Early Dark Energy, Phys. Rev. D 
84 (2011) 023518 [arXiv: 1104.0760 [astro-ph.CO]]. 

[12] E. Calabrese, M. Martinelli, S. Pandolfi, V. F. Cardone, C. J. A. P. Martins, S. Spiro and 

P. E. Vielzeuf, Dark Energy coupling with electromagnetism as seen from future low-medium 
redshift probes, Phys. Rev. D 89 (2014) 8, 083509 [arXiv:1311.5841 [astro-ph.CO]]. 

[13] L. Amendola, A. C. O. Leite, C. J. A. P. Martins, N. J. Nunes, P. O. J. Pedrosa and 

A. Seganti, Variation of fundamental parameters and dark energy. A principal component 
approach, Phys. Rev. D 86 (2012) 063515 [arXiv:1109.6793 [astro-ph.CO]]. 

[14] A. C. O. Leite, C. J. A. P. Martins, P. O. J. Pedrosa and N. J. Nunes, Fundamental Cosmology 
from Precision Spectroscopy: 1. Varying Couplings, Phys. Rev. D 90 (2014) 6, 063519 
[arXiv: 1409.3963 [astro-ph.CO]]. 


- 19 


[15] A. C. O. Leite and C. J. A. P. Martins, Phys. Rev. D 91 (2015) 103519 [arXiv:1505.05529 
[astro-ph.CO]]. 

[16] C. J. A. P. Martins and A. M. M. Pinho, Phys. Rev. D 91 (2015) 103501 [arXiv:1505.02196 
[astro-ph.CO]]. 

[17] G. R. Dvali and M. Zaldarriaga, Changing alpha with time: Implications for fifth force type 
experiments and quintessence, Phys. Rev. Lett. 88 (2002) 091303 [hep-ph/0108217]. 

[18] T. Chiba and K. Kohri, Quintessence cosmology and varying alpha. Prog. Theor. Phys. 107 
(2002) 631 [hep-ph/0111086]. 

[19] N. J. Nunes and J. E. Lidsey, Reconstructing the dark energy equation of state with varying 
alpha, Phys. Rev. D 69 (2004) 123511 [astro-ph/0310882]. 

[20] P. E. Vielzeuf and C. J. A. P. Martins, Varying constants and dark energy with the E-ELT, 
Mem. Soc. Ast. It. 85 (2014) 1, 155 [arXiv:1309.7771 [astro-ph.CO]]. 

[21] T. Damour and J. F. Donoghue, Phenomenology of the Equivalence Principle with Light 
Scalars, Class. Quant. Grav. 27 (2010) 202001 [arXiv:1007.2790 [gr-qc]]. 

[22] T. Damour, F. Piazza and G. Veneziano, Runaway dilaton and equivalence principle violations, 
Phys. Rev. Lett. 89 (2002) 081601 [gr-qc/0204094]. 

[23] C. J. A. P. Martins, P. E. Vielzeuf, M. Martinelli, E. Calabrese and S. Pandolfi, Evolution of 
the fine-structure constant in runaway dilaton models, Phys. Lett. B 743 (2015) 377 
[arXiv: 1503.05068 [astro-ph.CO]]. 

[24] H. B. Sandvik, J. D. Barrow and J. Magueijo, A simple cosmology with a varying Hne structure 
constant, Phys. Rev. Lett. 88 (2002) 031302 [astro-ph/0107512]. 

[25] P. M. M. Leal, C. J. A. P. Martins and L. B. Ventura, Fine-structure constant constraints on 
Bekenstein-type models, Phys. Rev. D 90 (2014) 2, 027305 [arXiv: 1407.4099 [astro-ph.CO]]. 

[26] N. Suzuki, D. Rubin, C. Lidman, G. Aldering, R. Amanullah, K. Barbary, L. F. Barrientos and 
J. Botyanszki et ai, The Hubble Space Telescope Cluster Supernova Survey: V. Improving the 
Dark Energy Constraints Above z^l and Building an Early-Type-Hosted Supernova Sample, 
Astrophys. J. 746 (2012) 85 [arXiv:1105.3470 [astro-ph.CO]]. 

[27] O. Farooq and B. Ratra, Hubble parameter measurement constraints on the cosmological 
deceleration-acceleration transition redshift, Astrophys. J. 766 (2013) L7 [arXiv: 1301.5243 
[astro-ph.CO]]. 

[28] T. Rosenband et ai, Frequency Ratio of A1+ and Hg+ Single-Ion Optical Clocks; Metrology at 
the 17th Decimal Place, Science, 319 (2008) 1808. 

[29] F. Luo, K. A. Olive and J. P. Uzan, Phys. Rev. D 84 (2011) 096004, [arXiv:1107.4154 [hep-ph]]. 

[30] M. C. Ferreira, M. D. Juliao, C. J. A. P. Martins and A. M. R. V. L. Monteiro, Phys. Rev. D 
86 (2012) 125025, [arXiv: 1212.4164 [hep-ph]]. 

[31] A. Songaila and L. L. Cowie, Constraining the Variation of the Fine Structure Constant with 
Observations of Narrow Quasar Absorption Lines, Astrophys. J. 793 (2014) 103 

[arXiv: 1406.3628 [astro-ph.CO]]. 

[32] P. Molaro, D. Reimers, 1. 1. Agafonova and S. A. Levshakov, Bounds on the fine structure 
constant variability from Fell absorption lines in QSO spectra, Eur. Phys. J. ST 163 (2008) 
173 [arXiv:0712.4380 [astro-ph]]. 

[33] H. Chand, R. Srianand, P. Petitjean, B. Aracil, R. Quasi and D. Reimers, On the variation of 
the fine-structure constant: very high resolution spectrum of qso he 0515-4414, Astron. 
Astrophys. 451 (2006) 45 [astro-ph/0601194]. 


- 20 - 


[34] 1. 1. Agafonova, P. Molaro, S. A. Levshakov and J. L. Hon, First measurement of Mg isotope 
abundances at high redshifts and accurate estimate of Delta alpha/alpha, Astron. Astrophys. 
529 (2011) A28 [arXiv:1102.2967 [astro-ph.CO]]. 

[35] T. A. Wagner, S. Schlamminger, J. H. Gundlach and E. G. Adelberger, Torsion-balance tests of 
the weak equivalence principle, Glass. Quant. Grav. 29 (2012) 184002 [arXiv:1207.2442 [gr-qc]]. 

[36] J. Muller, F. Hofmann and L. Biskupek, Testing various facets of the equivalence principle 
using lunar laser ranging, Glass. Quant. Grav. 29 (2012) 184006. 

[37] Z. Slepian, J. R. Gott, III and J. Zinn, Will Upcoming Observations Distinguish Between 
Slow-Roll Dark Energy and a Cosmological Constant?, Mon. Not. Roy. Astron. Soc. 438 
(2014) 1948 [arXiv:1301.4611 [astro-ph.CO]]. 

[38] D. J. E. Marsh, P. Bull, P. G. Ferreira and A. Pontzen, Quintessence in a quandary: Prior 
dependence in dark energy models, Phys. Rev. D 90 (2014) 10, 105023 [arXiv:1406.2301 
[astro-ph.CO]]. 

[39] J. D. Barrow, H. B. Sandvik and J. Magueijo, The behavior of varying alpha cosmologies, 

Phys. Rev. D 65 (2002) 063504 [astro-ph/0109414]. 

[40] J. Berg et al. [MICROSCOPE Collaboration], J. Phys. Conf. Ser. 610 (2015) 012009 
[arXiv:1501.01644 [gr-qc]]. 

[41] J. Overduin, F. Everitt, P. Worden and J. Mester, STEP and fundamental physics, Class. 
Quant. Grav. 29 (2012) 184012 [arXiv:1401.4784 [gr-qc]]. 

[42] E. D. Davis, C. R. Gould and E. I. Sharapov, Oklo reactors and implications for nuclear 
science, Int. J. Mod. Phys. E 23 (2014) 4, 1430007 [arXiv:1404.4948 [nucl-th]]. 

[43] Y. V. Petrov, A. I. Nazarov, M. S. Onegin, V. Y. Petrov and E. G. Sakhnovsky, Natural 
nuclear reactor Oklo and variation of fundamental constants. Part 1. Computation of neutronic 
of fresh core, Phys. Rev. C 74 (2006) 064610 [hep-ph/0506186]. 

[44] C. R. Gould, E. I. Sharapov and S. K. Lamoreaux, Time variability of alpha from realistic 
models of Oklo reactors, Phys. Rev. C 74 (2006) 024607 [nucl-ex/0701019]. 

[45] M. S. Onegin, Investigation of the fundamental constants stability based on the reactor Oklo 
burn-up analysis. Mod. Phys. Lett. A 27 (2012) 1250232 [arXiv:1010.6299 [nucl-th]]. 


- 21 - 


